Backscattering-Induced Crossover in Deterministic Diffusion 
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We study diffusion in a one-dimensional periodic array of scatterers modeled by a simple map. 
The chaotic scattering process for this map can be changed by a control parameter and exhibits the 
dynamics of a crisis in chaotic scattering. We show that the strong backscattering associated with 
the crisis mechanism induces a crossover which leads to different asymptotic laws for the parameter- 
dependent diffusion coefficient. These laws are obtained from exact diffusion coefficient results and 
are supported by simple random walk models. We argue that the main physical feature of the 
crossover should be present in many other dynamical systems with non-equilibrium transport. 
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One of the basic mechanisms in the theory of chaotic 
dynamical systems are so-called crisis events, where the 
asymptotic dynamics of the system changes drastically 
with respect to the variation of a control parameter 
10-^. Recently, it was discovered that similar events oc- 
cur in simple chaotic scattering systems when the scat- 
tering rules are varied. This phenomenon has been de- 
noted as a crisis in chaotic scattering On the other 
hand, over the past few years a considerable literature 
has been developed in which the origin of transport in 
non-equilibrium statistical mechanics has been related 
to the characteristics of chaotic scattering processes [||. 
One problem studied was deterministic diffusion in sim- 
ple one-dimensional maps |-§, where for some exam- 
ples parameter-dependent diffusion coefficients have been 
computed by taking the complete equations of motion of 
the dynamical systems into account PjIo|]. Similar one- 
dimensional maps have been proposed in Ref. Q as sim- 
ple models for a crisis in chaotic scattering. Thus, the 
aim of this Letter is to investigate whether a crisis in 
chaotic scattering has an impact on deterministic diffu- 
sion in simple one-dimensional maps. 

In the following, we consider discrete one-dimensional 
piecewise linear chaotic maps with uniform slope, x„+i = 
Mh{xn), where his a control parameter, and Xfi IS the po- 
sition of a point particle at the discrete time n. Mh (x) is 
continued periodically beyond the interval [0, 1] onto the 
real line by a lift of degree one, Mh{x -f 1) = Mh{x) -\- 1. 
We assume that Mh(x) is anti-symmetric with respect to 
X = 0, Mh(x) — —Mfi{—x), i.e., that there is no drift 
imposed on a point particle. As an example, we consider 
the sawtooth map sketched in Fig. |l|. It was chosen as a 
periodic continuation of the map studied in Ref. which 
exhibits a crisis in chaotic scattering. The control param- 
eter is here the height h of the map which is related to the 
absolute value of the slope ahy h — (a — 3)/4. The diffu- 
sive properties of similar maps have been studied in Refs. 
1^,^. For this sawtooth map the parameter-dependent 
diffusion coefhcient has been computed by solving the 
Frobenius-Perron equation of the dynamical system B] , 



Pn+i{x) ^ dy pn{y) S{x- Mh{y)) 



(1) 



where Pn{x) is the probability density for points on the 
real line, and Mh{y) is the map under consideration. 
There exists a dense set of parameter values h for which 
one can construct Markov partitions for the map, and for 
each of these parameter values Eq. (|l|) can be written as 
a matrix equation ||9|,[l0||. 



P„+i = (l/|a|)rp„ 



(2) 



p„ represents a column vector of the probability densi- 
ties in each part of the Markov partition at time n, and 
T is a topological transition matrix which can be ob- 
tained from the Markov partition. However, instead of 
solving the eigenvalue problem of T Q , here solutions for 
the probability density vector p„ have been obtained by 
iterating Eq. (||), 



Pn+i - (l/|ar)7^>o 



(3) 



Starting with any probability density vector pg this it- 
eration method enables us to compute the exact time- 
dependent probability density p„ at any time step n and 
all other dynamical quantities based on probability den- 
sity averages for maps of the type of Mh{x) pO[-p^. In 
particular, it provides an efhcient way to calculate diffu- 
sion coefhcients by employing an Einstein formula. 



D{h) = lim -!- f dx pn{a 

n^oo In J 



(4) 



where the integral is the second moment of the time- 
dependent probability density . 

Fig. |l| shows a log-log plot of the diffusion coefficient 
as a function of /i up to /i = 3.5. Included are four curves 
which describe the coarse-grained behavior of the exact 
results. For integer values of the height the diffusion 
coefficient can be computed analytically by applying the 
eigenvalue method of Ref. [p| , and we get 
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The two dashed curves give approximate limits for the 
oscillations of the exact diffusion coefficient in the range 
h > he- They are obtained by fitting the diffusion coeffi- 
cient with the functional form of Eq. at /i = (2fc + l)/2 
and h = {4k + 3)/4, G A^o, for the upper and lower 
curve, respectively. The two dotted curves show two sim- 
ple random walk approximations. For large heights the 
distance a point particle travels at one time step by mov- 
ing from one unit interval to another is taken into account 
exactly and we get 



Drwlih) 
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which gives the dotted line plotted for h > he- For small 
heights the absolute value of the distance is approximated 
to either zero or one, depending on whether the particle 
remains on a unit interval or moves from one unit interval 
to the next j?). This leads to 



Drw2{h) 



2h 



4/1 + 3 



0) 
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These approximations indicate three different regions of 
coarse-grained behavior for the exact diffusion coefficient: 
The first one is a simple initial region, where the diffu- 
sion coefficient behaves linearly for small heights. For 
h > he it decreases slightly on increasing the height. Fi- 
nally, for h > 0.5 it starts to grow quadratically in the 
height, but with strong oscillations on a fine scale. The 
transition between the two different types of asymptotic 
coarse-grained behavior, which occurs in the intermedi- 
ate region oi he < h < 0.5, can be understood by referring 



Q 0.1 



0.01 



0.01 




to the action of certain microscopic scattering mecha- 
nisms of the map, which are closely connected to the 
dynamics of a crisis in chaotic scattering. 

These scattering mechanisms are introduced in Fig. 
where certain regions of the map have been distinguished 
by shaded squares and triangles: The triangles refer to 
parts where points of one unit interval get mapped from 
that particular interval into another unit interval. Addi- 
tionally, if points enter a square they preferably move into 
the triangular escape region above or below the respective 
square after some iterations. These squares are identical 
to the squares of an analogous scattering model, where 
they provide the fundamental mechanism for a crisis in 
chaotic scattering [Q . The abbreviations / ( "forward" ) 
and b ( " backward" ) in these scattering regions refer to 
the dynamics of the critical point of the map, which is 
indicated by a small circle. Its first iteration is shown by 
the dashed line with the arrows. At its second iteration, 
and by increasing the height h continuously up from zero, 
the orbit of the critical point, denoted as the critical or- 
bit in the following, travels along the graph of the map 
in the next right box from the upper left to the lower 
right, as indicated by bold black arrows. This way, the 
critical orbit explores all the different scattering regions 
of the map as h is increased from zero. If it hits a region 
labeled by a 6 it is in a position to get backscattered into 
the box to the left. Vice versa, if the orbit enters an / re- 
gion it is in a preferable position to move further forward 
to the next box to the right. The critical point indicated 
in Fig. ^ is part of a forward scattering region. Note that 
there is a dense set of points around the critical point 
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FIG. 1. Double logarithmic plot of the diffusion coefficient 
D{h) with respect to the height h for the sawtooth map shown 
in the figure. The graph is based on 38,889 single data points. 
Two random walk solutions (dotted lines) and two curves 
which approximately give the boundaries of the osciUations of 
D(h) for values of h above the crisis point he (dashed hues) 
are included. 



FIG. 2. Chaotic scattering in the sawtooth map and its 
connection to the behavior of the diffusion coefficient D{h). 
Certain microscopic scattering mechanisms of the map are 
identified by shaded squares and triangles. The same sym- 
bols are shown along the D{h) curve, where they indicate the 
impact of the respective scattering regions on the diffusion 
coefficient. The graph consists of 10,268 single data points. 



2 



which exhibits the same dynamics, at least for the first 
few iterations. An event dynamically analogous to the 
one of a crisis in chaotic scattering now occurs at the 
crisis point, defined by the parameter value of the height 
he for which the critical orbit of the forward scattering 
region hits the boundary of a backward scattering region 
in the next right box after one iteration for the first time. 
This case is illustrated in Fig. ^ and, for the map shown, 
is determined by he — (-\/T7 — 3)/8 ~ 0.1404. We remark 
that this process is topologically not identical to a crisis 
in chaotic scattering in the sense that a crisis is generated 
by the merging of two formerly isolated invariant sets in 
the phase space. Nevertheless, we argue that dynami- 
cally this process provides the same characteristics as a 
crisis in chaotic scattering, especially the onset of strong 
backscattering 

The squares and the triangles along the diffusion coeffi- 
cient curve now refer to parameter regions where the crit- 
ical orbit gets mapped into the respective scattering re- 
gions. The different symbols on the curve denote bound- 
ary points where the critical orbit enters, or leaves, these 
regions. The diffusion coefficient clearly decreases glob- 
ally if the critical orbit enters a backscattering region, 
and it increases globally if it exhibits forward scatter- 
ing. Hence, the different microscopic scattering mech- 
anisms defined above are connected to regions in the 
macroscopic diffusion coefficient which exhibit different 
parameter-dependent behavior. Most importantly, the 
crisis point he corresponds to the first strong local max- 
imum of the curve, as shown in Fig. ^. Magnifications 
reveal that the fine structure of the curve changes signifi- 
cantly just below and above the crisis value, and that the 
curve is fractal This can be understood in detail by 
refining the procedure explained above and is reported 
elsewhere pO| , pl| . Moreover, in Fig. ^the crisis point he 
separates the regions which are described by the two dif- 
ferent random walk models and which show significantly 
different coarse-grained behavior. 

We conclude that the onset of strong backscattering 
affects the diffusion coefficient of the sawtooth map not 
only on a fine scale, but also on a coarse-grained scale. 
This phenomenon may be understood as a backscattering- 
induced crossover in deterministic diffusion. However, 
such a dynamical event must eventually take place in 
any map of the type of M}^[x) at a certain parameter 
value, independently of its special functional form and 
independently of the special topological characteristics 
of a crisis in chaotic scattering. This can be checked by 
identifying the forward and backward scattering regions 
of a map and applying the definition of the crisis point 
given above. As an example, the piecewise linear, discon- 
tinuous, non-sawtooth map studied in Ref. |^ has been 
analyzed. It turns out that the respective crisis point of 
the map again corresponds to the first strong local maxi- 
mum of the diffusion coefficient curve, and that again the 
crisis point is related to a change between two different 
laws for the parameter-dependent diffusion coefficient on 
a coarse-grained scale: As before, the diffusion coefficient 



grows linearly for small values of the height, although 
with a slope different from that obtained for the saw- 
tooth map. It lacks a broad crossover region right above 
the crisis point, but as in case of the sawtooth map it in- 
creases quadratically with a factor of 1/6 in the limit of 
large heights. Analogous results for asymptotic diffusion 
coefficients have been reported in Ref. for a variety 
of other maps of the type of Mfi{x), although in this 
previous work diffusion coefficients could be computed 
exactly only for special values of the height. This leads 
us to conjecture that a backscattering-induced crossover, 
as described above, is typical for diffusive maps of the 
type of M/j(x). Moreover, we suspect that for these 
maps the parameter-dependent coarse-grained diffusion 
coefficient always decreases linearly in the limit of small 
heights, and that it increases quadratically in the limit of 
large heights with a universal factor of 1/6. Similar re- 
sults concerning the asymptotic coarse-grained behavior 
of diffusion coefhcients may be found in certain classes of 
two-dimensional maps, as, e.g., in standard and sawtooth 
maps, where quadratic laws in the limit of large control 
parameters have already been obtained P,p^. 

We suppose that the main physical feature of the 
crossover, i.e., the connection between the onset of strong 
microscopic backscattering and a change in the behav- 
ior of macroscopic parameter-dependent transport coef- 
ficients, is quite common not only in discrete one- and 
two-dimensional models, but also in higher-dimensional 
Hamiltonian systems: For example, in Ref. ||l^ diffusion 
of an electron in a two-dimensional crystal, modeled by 
Coulombic periodic potentials, has been studied. For this 
class of dynamical systems an energy threshold has been 
proved to exist above which the diffusion coefhcient in- 
creases with a power law in the particle energy, whereas 
below this threshold no diffusion coefficient exists. In 
Ref. Q it has been found that related models exhibit a 
crisis in chaotic scattering. Hence, we conjecture that 
the existence of this energy threshold is linked to the dy- 
namics of a crisis in chaotic scattering as in case of the 
backscattering-induced crossover discussed above. Fur- 
thermore, for the models with a crisis in chaotic scatter- 
ing the significance of orbiting collisions has been pointed 
out, which indicate the onset of strong backscattering [||. 
However, orbiting collisions have already been discussed 
extensively for fluid systems at low densities, where parti- 
cles interact via Lennard- Jones potentials, and a qualita- 
tive connection between the onset of these collisions and 
a small change in the temperature-dependent behavior 
of transport coefRcients has been established Thus, 
physically a direct line may be drawn from the dynamical 
origin of crises in dynamical systems over the occurrence 
of certain microscopic chaotic scattering processes to a 
specific coarse-grained behavior of transport coefRcients, 
which may be related to macroscopic dynamical crossover 
phenomena, or in certain cases possibly even to dynami- 
cal phase transitions. 

We conclude with a few remarks: (a) The sawtooth 
map under consideration here was chosen as a diffusive 
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version of a one-dimensional dynamical system which ex- 
hibits a crisis in chaotic scattering. A crossover in the 
parameter-dependent diffusion coefficient of this map has 
been found which is generated by the dynamical mecha- 
nism of this crisis event. It affects the diffusion coefficient 
not only on a fine scale, but also on a coarse-grained scale 
and leads to two different algebraic laws for the asymp- 
totic diffusion coefficient, (b) These different algebraic 
laws, i.e., a linear decrease of the diffusion coefficient for 
small heights and a quadratic increase with a factor of 
1/6 in the limit of large height, are suspected to be uni- 
versal for maps of the type of Mh{x). (c) The crossover 
is understood physically by relating the onset of strong 
microscopic backscattering to a change in the behavior 
of the macroscopic parameter-dependent transport coef- 
ficient. This main physical feature is conjectured to be 
quite common in dynamical systems which exhibit non- 
equilibrium transport. 
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